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Abstract 

We examine five years of rain rate and modeled slant path attenuation distributions at 
20 GHz and 30 GHz derived from a network of 10 tipping bucket rain gauges. The rain 
gauge network is located within a grid 70 km north-south and 47 km east-west in the Mid- 
Atlantic coast of the United States in the vicinity of Wallops Island, Virginia. Distributions 
were derived from the variable integration time data and from one minute averages. It 
was demonstrated that for realistic fade margins, the variable integration time results are 
adequate to estimate slant path attenuations at frequencies above 20 GHz using models 
which require one minute averages. An accurate empirical formula was developed to convert 
the variable integration time rain rates to one minute averages. Fade distributions at 20 GHz 
and 30 GHz were derived employing Crane’s Global model because it was demonstrated to 
exhibit excellent accuracy with measured COMSTAR fades at 28.56 GHz. 

1. Introduction 

Earth-satellite commur ications at frequencies above 10 GHz suffer from attenuation 
caused by rain. Designers of such systems are interested in having a priori knowledge of 
the probability of exceeding different levels of rain attenuations so as to establish appropri- 
ate fade margins into their systems. They are also interested in establishing estimates of the 
year to year variability of rain fade margins for particular geographic regions so that commu- 
nication systems reflect the extremes of these variabilities. Direct measurements of beacon 
signals from geostationary satellites have been a means to determine the above information 
fCCIR 1986; Goldhirsh, 1982] and experiments are presently being pursued in Europe an 
the United States with satellites such as Intelsat [Vogel and Torrence, 1991] and Olympus 
[Satellite Communications Group, 1991], and in the near future with ACTS [Davarian, 1991]. 
Such measurements are also important for model development and validation. 

The employment of rain-fade models is an ancillary approach for arriving at the above 
information. An important class of models use rain rate data acquired from rain gauges in 
different geographic regions for estimating cumulative distributions of slant path attenuation 
[CCIR 1986]. Implementation of these models has the advantage that is relatively inexpen- 
sive, may be made over many years, and may address questions not easily addressed with 
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direct measurements; sucli as “What is the variability of fade distributions over distances of 
10s of kilometers within a given geographic region?” and “What is the network average year 
to year variability in fades?” It is the objective of this paper to address the above questions 
through the presentation of analytical results of five years of ram rate data from a network 
of 10 rain gauges in the Mid- Atlantic coast of the United States. The results presented here 
is an expansion and elaboration of a previous effort by Goldhirsh (1990], encompassing two 
years of measurements. 

Many slant path attenuation models employ one minute averages in their models [CCIR, 
1986]. This creates a complication for investigators who utilize “tipping” bucket rain gauges 
because the rain rates are measured with variable integrations times. It is therefore another 
objective of this effort to assess the sensitivity to estimation errors of slant path attenuation 
distributions derived by variable integration times obtained from tipping bucket rain gauges. 
This question will be examined by comparing the variable integration time results with those 
derived employing one minute rain rate averages. 

2.0 Description of Rain Gauge System 

The rain gauges have an eight inch diameter collecting cylinder and contain buckets 
which tip after an accumulation of 0.254 mm of rainfall. After each tip, a switch at the 
gauge closes for approximately 100 milliseconds (switch closure time). The switch closure 
changes a voltage level monitored by a connecting PC. Whenever this voltage level change 
is noted, the PC records the computer clock time. Every two hours the accumulated tipping 
times are automatically recorded on files on a 5 1/4” disk. Also stored on the disk are the 
Julian date, the local time the file was written onto disk, and the total rainfall in mm. 

The integration time AT (time between tips) is variable for the tipping bucket rain gauge 
and is defined by 

R = (mm/h) l 1 ) 

AT v 

where AT is expressed in seconds (integration time) and R is the rain rate in mm/h. We 
note that a 60 second integration time corresponds to a rain rate of 15.2 mm/h. Hence, 
measured rain rates smaller than approximately 15 minutes have integration times which 
are longer than one minute. A rain rate of 1 mm/h has an integration time of 15 minutes 
which represents approximately the lower measurement threshold with the above described 

gauge. 

We report on the results of 10 rain gauges located within a griddcd region of 70 km 
north-south and 47 km east-west in the Mid-Atlantic coast of the United States. All the 
gauges are located within a radial distance of 60 km from the SPANDAR radar facility at 
the NASA Goddard Wallops Flight Facility (WFF), Wallops Island, Virginia. The 10 gauge 
locations are depicted in the map of Figure 1. The individual rain gauge systems are located 
at the home grounds of staff working at the WFF. The staff maintains these systems on a 
continuous basis. The floppy disks are removed on a weekly basis and submitted to a central 
processor for reduction and analysis. Careful calibrations for each of the site gauges are 
performed twice per year and the system is maintained with errors of less than 5% in ram 
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rate at rates of 12 to 15 mm/h [Gebo, 1991]. 

3.0 Network Rain Rate Distributions 

The rain rate distributions described in the following paragraphs were derived from mea- 
surements of 10 rain gauge sites encompassing the period June 1, 1986 through May 31, 
1991 with the following caveats; the data related to Site #4 were only available for four 
years (June 1, 1987 through May 31, 1991), and the data for Site #6 for three years (June 1, 
1986 — May 31, 1989). The time and network average distributions were weighted to reflect 
the shorter operational time periods of these two sites. 

3.1 Overall Average Rain Rate Distribution 

In Figure 2 is plotted the combined distribution comprising the spatial (10 sites) and 
temporal (five years) average rain rate distributions for the variable integration time. Also 
plotted for comparison is the one minute integration time case. The combined average will 
hereafter be referred to as the “overall average” case. The one minute distribution was 
obtained by averaging the variable rain rate-time series over one minute contiguous periods 
and determining the cumulative distributions for the one minute average rain rate-time 
series. These two rain rate distributions are noted to be different by less than 3 mm/h up 
to 60 mm/h (0.01 % probability) where the time between tips is approximately 15 seconds. 

It was shown by the authors [Goldhirsh et ah, 1992] that at the 0.05% level, the differences 
between the modeled attenuations derived from the variable integration time and the one 
minute average distributions are negligibly small at both 30 GHz (0.6 dB) and 20 GHz (0.3 
dB) compared to the respective fade levels (31.1 dB and 13.8 dB). We therefore conclude 
that it is not required to use one minute averages in the modeled rain rates over the range 
of practical (modeled) fade margins at 20 GHz and 30 GHz. 

3.2 Scaling Formulation Rain Rate Distributions 

For completeness, however, we address the question as to how one may scale the rain 
rate distributions obtained with the variable integration time gauge to distributions which 
correspond to one minute averages. This formulation may be used to assess attenuations at 
lower frequencies where higher rain rates may define realistic fade margins. Thereafter, rain 
rate and modeled attenuation distributions will be presented here which correspond only to 
the variable integration time rain rate results. 

A formulation has been derived which converts the variable rain rate distributions to the 
one minute distributions in the percentage interval 0.1% to 0.001% [Goldhirsh et ah, 1992]. 
At percentages greater than 0.1%, the distributions are virtually identical for the described 
rain gauge. The formulation is given by 

R!(P) = R v (P)-6(P) (2) 

where 

6(P) = a + 0 R + 7 R 2 (3) 
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( 4 ) 


( a = - 0.2036 
l f3 = 2.250 x 10~ 2 
( 7 = 4.729 x 10" 4 

and where Ri(P) is the one minute averaged rain rate (mm/h), R«(P) is the variable inte- 
gration time rain rate, and ^(P) is the correction factor given by (3) and (4). All of these 
quantities are taken at the probability, P. The above formulation predicts the 60 second 
“overall average” distribution to within 2.5 mm/h at the 0.001% level. 

3.3 Temporal Variability-Rain Rate Distributions 

To establish a measure of the year to year variability in the rain rate distributions, 
we show in Figure 3 five distributions corresponding to the network average for each of 
five years (1986-87, 87-88, 88-89, 89-90, 90-91; hereafter referred to as years 1 through 5, 
respectively). Also plotted (dotted curve) is the “overall average” distribution. We note 
from Figure 3 that four of the years for the network average show similar distributions and 
year 4 shows considerably larger rain rates. For example, at 0.1%, year 4 shows a rain rate 
of approximately 24 mm/h compared to a range of other rain rates between approximately 
11 mm/h and 15 mm/b. This result is a demonstration of the need to obtain multi-year 
rain rate distributions in order to assess extreme levels that may arise. At the 0.1% level a 
maximum rain rate difference of 13 mm/h exists. This is the difference between values for 
year #4 (highest rain rate) and year #5 (lowest rain rate). 

3.4 Spatial Variability-Rain Rate Distributions 

In this section we examine the variability in the distributions caused by measurements at 
different locations separated by 10s of kilometers within the Mid- Atlantic coast geographic 
region. In Figure 4 is given the rain rate distributions for each of 10 sites averaged over 
the five year period; with the exception of Sites #6 and #4 averaged over 3 and 4 years, 
respectively. It is apparent that the variability in the distributions for the individual site 
locations (Figure 4) appears to be smaller than the year to year variability (Figure 3). The 
maximum rain rate differences at 0.1%, 0.01%, and 0.001% are 3.5, 12.2, and 18.0 mm/h, 
respectively. Shown also, (dotted curve) is the overall network and temporal average. 

We observe from Figure 4 that Site $1 shows generally larger values, whereas Site $6 
shows generally minimal levels of rain rates. Site #6 (Figure 1) is located closest to the 
shoreline, where the rains tend to diminish because of the cooling effects of the sea water. 
That is, the cooler sea water tends to dampen the more intense convective rain cells which 
originate over land and are sustained by ground heating, as for example, air mass type 
systems. Site #1 is the most northerly site located approximately 10 km from the shore but 
furthest from the Chesapeake Bay than the other sites. Storm systems, which are generally 
out of the south-west, w ill move over a greater land mass and be sustained by ground heating 
as they pass Site #1 as compared to the other sites. 
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4.0 Attenuation Distributions 

4.1 Employment of the Global Model 

We use here the Crane’s Global model [Crane, 1980] to derive the corresponding attenu- 
ation distributions at 20 GHz and 30 GHz at a slant path elevation angle of 45°. This model 
was employed here because of it’s success in estimating the COMSTAR, derived attenuation 
distribution at 28.56 GHz [Goldhirsh, 1982] when applying the measured rain rate distribu- 
tion [Goldhirsh, 1990; Goldhirsh et al, 1992]. A comparison between the directly measured 
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RAIN RATE (inm/h) 

Figure 3: Comparison of yearly network average rain rate distributions over years 1 through 

5. 



RAIN RATE (rnm/h) 

Figure 4: Comparison of rain rate distributions for 10 sites. Each site was averaged over five 
years with the exception of Site #6 (3 years) and Site #4 (4 years). 
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Figure 5: Comparison of yearly network average attenuation distributions at 20 GHz over 
years 1 through 5. 



ATTENUATION (dB) 

Figure 6: Comparison of yearly network average attenuation distributions at 30 GHz over 
years 1 through 5. 
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e 7: Comparison of attenuation distributions at 20 GHz for 10 sib 
ged over five years with the exception of Site #6 (3 years) and Site 
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Figure 8: Comparison of attenuation distributions at 30 GHz for 10 sites. Each site was 
averaged over five years with the exception of Site #6 (3 years) and Site #4 (4 years). 


